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ABSTRACT 

We study the interplay between gas phase metalHcity (Z), specific star formation rate 
(SSFR) and neutral hydrogen gas (HI) for galaxies of diff'erent stellar masses. Our study uses 
spectroscopic data from GAM A and SDSS star forming galaxies, as well as Hl-detection 
from the ALFALFA and GASS public catalogues. We present a model based on the Z-SSFR 
relation that shows that at a given stellar mass, depending on the amount of gas, galaxies will 
follow opposite behaviours. Low-mass galaxies with a large amount of gas will show high 
SSFR and low metallicities, while low-mass galaxies with small amounts of gas will show 
lower SSFR and high metallicities. In contrast, massive galaxies with a large amount of gas 
will show moderate SSFR and high metallicities, while massive galaxies with small amounts 
of gas will show low SSFR and low metallicities. Using ALFALFA and GASS counterparts, 
we find that the amount of gas is related to those drastic diff'erences in Z and SSFR for galaxies 
of a similar stellar mass. 

Key words: galaxies: abundances, galaxies: fundamental parameters, galaxies: star forma- 
tion, galaxies: statistics 



1 INTRODUCTION 

The formation of galaxies is intimately dependent on the conver- 
sion of stars into gas, the production of heavy elements, recycling 
of this material into the interstellar medium, and repetitions of this 
cycle. A detailed understanding of the interplay between each of 
gas mass, star formation rate, and metallicity is clearly important to 
understand the galaxy evolution process. Scaling relations between 
SFR and stellar or gas mass, and between mass and metallicity, 
have been explored for many years. More recently, the connection 
between gas fraction and metallicity have also begun to be explored 
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[Hughes et al.' (zOlZ). Identifying the causal relationships between 
these empirical scaling relations is, however, more challenging. 

From the fundamental properties of galaxies, we can identify 
extensive (scale-dependent) and intensive (scale-invariant) prop- 
erties. Intensive properties are essential or inherent to a system 
and are independent of the amount of baryonic mass, for example 
metallicity, temperature, density, and SSFR. Extensive properties 
on the other hand, depend on the amount of baryonic mass of the 
system, such as the star formation rate (SFR) and stellar mass (M^). 
By way of illustration of this point, if we were to split a galaxy in 
half, the metallicity of the halves are unchanged from the original 
galaxy. The SFR and stellar mass of the halves, however, are half 
that of the original, corresponding to the amount of baryonic mass 
removed. 
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Although the metalUcity is scale invariant, or intensive, it is 
relat ed to the stellar mass through the well known M^-Z relation 
fe.g. lTremonti et al. "20041: iLara-Lopez et alJl2009allbl) . In this rela- 
tionship, the metallicity of high mass systems is higher than that 
of low mass systems. This does not mean that the metallicity and 
M^ are dependent variables. A correlation or trend between an in- 
tensive and an extensive property cannot be directly causal. For ex- 
ample, simply combining a large number of low-metalliticy dwarf 
galaxies will increase the total mass, but will not directly increase 
the metallicity of the total system. Any such relation, consequently, 
must imply a common underlying cause, in this example likely to 
be the star formation history of the galaxy. 

On the other hand, relationships such as Z-SFR & Z-SSFR are 

I II 1 

shown to have w eaker co rrelation s (e.g. Lara-Lopez et al. 2010a; 

Lopez- SanchezI |201Q| : iLara-L opez, Lopez- Sanchez & Hopkins 



2 SAMPLE SELECTION 



2.1 Optical data 



2013h . Specifically, the Z-SSFR relation is formed by intensive 



properties of galaxies, and has not been widely studied since it 
shows high scatter. Nevertheless, metallicity is important for the 
ability of gas clouds to form stars, since it enables most of the gas 
cooling that will facilitate the star formation (Lebouteiller et al.l 
l2013h . 

Crucial information regarding galaxy properties is given 
by their neutral gas content. Large HI surveys such as the 
Arecibo Legacy F ast Arecibo L-band Feed Array (ALFALFA, 



Havnes et al. 2011), and the Galex Arecibo SDSS Survey (GASS, 
Catinella et a l. 2010), have provided new scaling relationships and 



dependencies for the cool gas in galaxies. 

The link bet ween the gas metalli city and HI content in galax- 
ies was studied in lzhang et alJ (120091) . Specifically, they study the 
HI dependence on the M^-Z relation using SDSS galaxies. They 
estimate empirically the HI content for 10^ emission-line galaxies 
in the SDSS-DR4, and fi nd that gas-poor gal axies are more metal 
rich at fixed stellar mass. iHughes et al.l (120121) confirms this result, 
from an investigation of the role of cold gas and environment on the 
M^-Z relation for 260 nearby late-type galaxies. They find that, at 
fixed stellar mass, galaxies with lower gas fractions typically also 
possess higher metallicities. In general, they observe that gas-poor 
galaxies are typically more metal rich, and demonstrate that the re- 
moval of gas from the outskirts of spirals increases the observed 
average metallicity by ~ 0. 1 dex. 

Athough some scaling relations with HI content have been 
studied, the interplay between metallicity and SSFR has not yet 
been analyzed as a function of the gas content or stellar mass. The 
Z-SSFR relation is formed by intensive properties, and is a key re- 
lationship to study the diff'erent properties of high- and low-mass 
galaxies. 

The gas mass is a direct measure of the available fuel in galax- 
ies to form stars. The relation of this gas with metallicity, stellar 
mass, and SSFR can tell us how fast a galaxy assembled stars in 
the past, the amount of stars it formed, and how actively it is form- 
ing stars at present. The main goal of this paper is to analyse the 
interplay between all these properties to produce a general picture 
of the gas recycling process. 

This paper is organized as follows. In §[2] we detail the data 
used for this study, and in §[3] we analyse the Z-SSFR relation 
and present a cartoon model based on it. In §|4l we study scaling 
relations using HI data. Finally, in §[5] we present our discussion 
and conclusions. Throughout we assume Hq - 70kms~^ Mpc~\ 
Qm = 0.3, Qa = 0.7. 



We consider data for emission-line galaxies from two large sur- 
veys, the Galaxy and Mass Assembly (GAM A phase-I survey, 
iDriver et al.l201 1^. and the Sloan Digital Sky Survey-Data Release 
7 (SDSS-DR7,[Abazajian et al. 2009). 

Data from the SDSS were taken with the 2.5 m telescope 
located at Apache Point Observatory (iGunn et aLll2006h . We use 
the emission-line analysis of SDSS-DR7 galaxy spectra performed 
by the MPA-JHU databas^j. From the full dataset, we only con- 
sider objects classi fied as galaxies in the "main galaxy sample" 
JStrauss et al.ll2002h with apparent Petrosian r-magnitudes in the 
range 14.5 < m^ < 17.7 7 and z < 0.33 . We u se gas metallicities 
measured as described in Trenionti et aLl(l2004h . SFR estimates de- 
scribed in Brinchmann et al. (2004), and total stellar masses esti- 
mated as in Kauff'mann et al. ( 2003a) . 

The GAMA phase-I is a spectroscopic survey with the 
3.9m Anglo- Australian Telescope (AAT) using the 2dF fibre 
feed and AAOmega multi-object spectrograph. For full de- 
tails o f the survey selection and properties see [Driver et al] 
(l201lh . SFRs measu rements are based on the Ha emission 
line as described in iGunawardhana et al.l (1201 lb . Metallicities 
were estimated using the empirical calibration provided by 
IPettini & Pagell (2004) between the oxygen abundance and the 
03N2 ( ([Oiii] ^5007/H/?)/([Nii] i6583/HQf) ) index. Both SFRs 
and metallicities were recalibrated to the Bayesian sy s tem de - 
scribed above using the calibrations of iLara-Lopez et al. (l2013bl). 
Finally, stellar masses were measured as described in iTavlor et al.l 
(2011). 

For both surveys, we selecte d only SF galaxies using the 
standard spectroscopic diagnostic (IBaldwin. Phillips & TerlevichI 
Il98lh . and using the discrimination of iKauff'mann et al.l (l2003bh . 
For reliable metallicity and SFR estimates, we selected galaxies 
with a signal-to-noise ratio (SNR) of 3 in Hor, Yip, and [N ii]. Galax- 
ies were selected in volume-limited samples in redshift and Pet- 
rosian r -band absolute magnitude as described in lLara-Lopez et al] 
(l2013bh . Our final SF optical sample is 35212 galaxies for GAMA 
and 156910 for SDSS. These samples are used to match with AL- 
FALFA and GASS data as described in the next section. 



2.2 HI data 

To test the correlations between optically deriv ed properties and the 
HI content of galaxies , we used the A LFALFA (IHavnes et alj201lh 
and the GASS ( Catinella et al.ll201Qh Surveys. 

ALFALFA is a blind survey of 21 cm HI emission over ~ 2800 
deg^ of sky. The public cataloguqj provides ~ 15,855 HI detec- 
tions, of which ~ 15,041 are associated with extragalactic objects. 
Since ALFALFA is a blind survey, we select only galaxies with 
Code=l, which refers to robust, reliably detected sources. We also 
remove sources with heliocentric velocities V/je//o< 100.0, which are 
unlikely to be galaxies. From this subsample, we cross-match RA, 
DEC, and redshift for our optical SF-sample described in § 12.11 ob- 
taining a final sample of 4443 SF galaxies with ALFALFA counter- 
parts. These come entirely from the SDSS, since the GAMA survey 
regions are not well covered by the publicly available ALFALFA 



^ http://www.mpa-garching.mpg.de/SDSS/ 

^ http ://egg . astro . comell . edu/alfalf a/data/index .php 
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data. Together with the fainter magnitude range of the GAM A tar- 
gets, this biases against them being identified with ALFALFA de- 
tections. 

In order to increase our HI coverage of massive galaxies, we 
match our SDSS SF-sample described in § l2.1l with the GASS pub- 
Uc catalogue. The GASS survey aims to observe a sample of ~ 
1000 galaxies selected from the SDSS spectroscopic and GALEX 
imaging surveys. GASS galaxies are selected to have stellar masses 
greater than W^ M^ and redshifts 0.025< z <0.05. The pubUc cat- 
alogue includes 232 galaxies, from which we find 48 counterparts 
with our optical SF-sample described in § 12.11 The rest of the 
GASS galaxies correspond to composite and AGN galaxies. The 
stellar mass range of our final sample is 10^ <M^< 10^^ Mq, with 
a median of ~ lO'^^M©. 



3 THE Z-SSFR RELATION 

The correlation between the amount of metals and the SSFR has 
not been explored in detail. Since both Z and SSFR are intensive 
properties, this relationship provides a direct way to analyse the 
diff'erent properties of low- and high-mass galaxies. 

Here we explore the median SSFR in bins of Z for diff'er- 
ent M^. In Fig. [T] we show the Z-SSFR relation for SDSS and 
GAMA data. Each of the panels corresp ond to diff'erent redshifts o f 
volume-limited samples as described in lLara-Lopez et al.1(l20l3b! ). 
The coloured circles correspond to the median SSFR in bins of Z 
for diff'erent M^. Opposing trends are clearly observed in Figs. [TJ 
and[TJ) between low- and high-mass galaxies with an inflection 
point at log(SSFR) ~ -9.9 yr~^ The SSFR for low-mass galax- 
ies shows an anticorrelation with Z while the SSFR for high-mass 
galaxies correlates with Z. 



3.1 A Model 

We suggest that the opposite behaviors for low and high mass 
galaxies observed in Fig.[T]can be understood assuming that galax- 
ies at a given M^ have diff'erent amounts of HI. Based on that figure, 
we generate the cartoon model shown in Fig.|2] In this model, low- 
mass galaxies are represented by a blue ellipse that extends from 
high SSFR and low Z, to low SSFR and high Z. Massive galaxies, 
represented by a red ellipse, extend from low SSFR and low Z to 
high SSFR and high Z. 

Diff'erent physical explanations can be given to justify this ba- 
sic supposition of our cartoon model. The presence of galactic in- 
flows and outflows and how they aff'ect the metal enrichment in 
galaxies is still a matter of debate. Galactic outflows have been 
shown to be importan t to reprod uce the M^-Z relation of galaxies 
using analytical (e.g., JErbluOOSh and hydro dynamic models ( e.g., 
iFinlator & Davil2008h . On the other hand, .Calura et all (l2009h re- 
produced the M^-Z relation mainly by means of an increasing ef- 
ficiency of star formation, without any need to invoke galactic out- 
flows. 

A general picture of downsizing suggests that low-mass 
galaxies process their gas slower than massive galaxies, therefore 
low-mass galaxies are actively forming stars, and present high SS- 
FRs today compared to massive galaxies. Hereafter, we will focus 
on downsizing to provide a possible explanation for the trends be- 
tween the intensive properties of SSFR and Z as a function of the 
extensive property of mass. 

A combination of downsizing and the amount of HI can ex- 
plain the opposing trends for low and high mass galaxies observed 



in the Z-SSFR relation. Focusing our attention on the low-mass 
galaxy branch (blue ellipse of Fig. |2]), we can explain the diff'er- 
ences in SSFR through the amount of HI present in each galaxy. 
Low-mass galaxies with a high amount of neutral gas will show 
higher SSFR today, as they have more fuel for star formation, than 
galaxies of the same mass with a lower amount of neutral gas. The 
same is applicable to high-mass galaxies (red ellipse of Fig. [21), in 
this branch, galaxies with a high amount of HI show higher SSFRs 
than galaxies at the same mass with lower amount of HI. 

It is noteworthy however, that Z plays the opposite role. While 
low-mass galaxies with a large amount of HI show low Z, massive 
galaxies with a large amount of HI show high Z. This can again be 
explained by downsizing, but not downsizing in stellar mass, rather 
the amount of HI is driving the rate of metal enrichment in galaxies. 
In the low-mass branch, galaxies with large HI show low Z because 
they are processing their gas slower and on longer time scales than 
galaxies with lower HI for the same stellar mass. On the other hand, 
massive galaxies with large HI show higher metallicities because 
they processed their gas faster in the past and have already reached 
high Z, and due to their large HI, they also have a high SSFR too. 



4 HI SCALING RELATIONS 

To test our model we use the 4491 SDSS optical counterparts with 
direct HI mass measurements from the ALFALFA and GASS sur- 
veys described in § 12.21 For the relationships shown in this section, 
we define the gas mass as Mgas = 1.32xMhi, and the gas mass 
fraction as Mgas / (M^ + Mgas). 

We generate the Z-SSFR relation and show it as a function of 
the gas mass fraction in Fig. |2] (right). Although with our HI sample 
it is only possible to see the low-mass branch in this relationship, it 
is clear that the gas mass fraction of galaxies increases as the SSFR 
increases, consistent with our proposed cartoon model. Black trian- 
gles in Fig. [21 (right) correspond to SDSS galaxies in the ALFALFA 
and GASS fields with no HI detection and M^>10^^ M©. A linear 
fit to these data is given by log(SSFR) = -3.509 - 0.7184 x, where 
x=12+log(0/H). 

Fig. [21 also suggests that there is a relation between the gas 
mass fraction and the SSFR. This indicates, regardless of the M^, 
that the intrinsic rate at which a galaxy is forming stars depends 
strongly on the gas mass fraction. This relationship is shown in Fig. 
[3^. Galaxies with a high gas mass fraction show a very steep rela- 
tion with SSFR. On the other hand, galaxies with a low gas mass 
fraction have uniformly flat and low SSFRs in which the lack of 
gas would prevent galaxies from forming stars. Indeed, there is an 

inflection point between this abrupt change at log(SSFR) 9.9, 

in agreement with the dividing line seen in Fig. [T] that motivated 
the cartoon model of Fig. [21 

The relations between the gas mass fraction with the stellar 
mass and metallicity are shown in Figs. [3]d and [3j:. Similarly to 
Fig.[3t, these relations show steep and flat tendencies for high and 
low gas mass fractions, respectively. Galaxies with high gas mass 
fractions have low metallicities and low-masses. This suggests, as 
downsizing indicates, that low-mass galaxies assemble their stars 
on longer timescales compared to massive galaxies. Galaxies with 
very low gas mass fractions, have assembled their stellar mass, ex- 
hausted their gas, and reached high stellar masses and high metal- 
licities. As a result, galaxies with low gas mass fraction will show a 
flat relation against stellar mass and metallicity, as seen in Figs.[3j) 
and[3l:, respectively. To further refine this picture, HI measurements 
for massive galaxies are necessary. 
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Figure 1. Z-SSFR relation for SDSS and GAMA data.The colour coded circles correspond to the median SSFR in bins of Z for different bins of M^. The grey 
points correspond to the density of SDSS and GAMA samples. 



CO 
CO 



decreasing SF efficiency 



moderate 
gas mass fraction 




i 



low gas 
I mass fraction 



moderate 
gas mass fraction 



Increasing SF efficiency 






CD 

cn 



O 



-lOr 



-11 7 



-0.8 -0.6 -0.5 -0.3 -0.2 0.0 






12+log(0/H) 



8.0 



8.5 
12+log(0/H) 



9.0 



Figure 2. Left: cartoon model of the Z-SSFR relationship based on the extreme stellar mass cases of Fig. \T\ The model is based on local redshift galaxies, 
and represents a snapshot of the behaviour of galaxies today. The blue and red ellipses represent low- and high-mass galaxies, respectively. Right: Z-SSFR 
relation for galaxies with ALFALFA and GASS counterparts. Galaxies are colour coded from low (magenta) to high (green) gas mass fraction. Black triangles 
show the median data for SDSS galaxies with no HI detection in the ALFALFA and GASS fields for log(M^) > 1 1.0 dex. 
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Figure 3. Scaling relations for SDSS galaxies with ALFALFA and GASS counterparts, (a) The gas mass fraction and SSFR relation colour coded from low 
(green) to high (magenta) metallicity. (b) The gas mass fraction and M^ relation colour coded as in Figure (a), (c) The gas mass fraction and Z relation for the 
same sample colour coded from low (magenta) to high (green) SSFR. 
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Another way to see the influence of the gas fraction on metals 
is through the M^-Z relation shown in Fig. [4^. Since low-mass 
galaxies are actively forming stars they show a high HI gas fraction 
but low metallicity, while massive galaxies have finished forming 
stars, exhausting their HI, and reaching high metallicities. A fit to 
this relation gives: 12+log(0/H) = 47.751 - 13.903 x+ 1.5839 x^ - 
0.05815 x\ with x=log(M^/Mo). 

Two important extensive variables analyzed here are the stellar 
mass and the total HI in galaxies. The relationship between both is 
shown in Fig. |41d. It is interesting to note that low-mass galaxies 
have a larger amount of mass in HI than in stars, which supports 
the scenario described in this section. Galaxies with a higher mass 
in stars were more efficient at forming stars in the past, resulting in 
a current amount of HI that is relatively smaller, although still larger 
in an absolute sense, than in low- stellar-mass systems. A fit to this 
relation gives: log(MHi/Mo) = -8.368 + 3.4166 x-0.1588 x^, with 
x=log(M,/Mo). 



5 DISCUSSION AND CONCLUSIONS 

A key question to answer to have a general picture of this model 
is, why do low metallicity galaxies have either high SSFRs or low 
SSFRs but almost nothing in between?. In the picture presented 
here, galaxies evolve in a downsizing fashion with two independent 
mechanisms that lead to having a galaxy with low metallicities: (/) 
A high-mass galaxy will be at low metallicity today if it exhausted 
its available fuel for star formation and can no longer continue to 
enrich its ISM. (//) A low-mass galaxy will be at low metallicity 
while it slowly forms stars, taking longer to enrich its ISM. 

These different scenarios lead naturally to the dearth of mid- 
range SSFRs for the very low metallicity systems. Staying at low-Z 
for massive galaxies means less and less star formation the lower 
the metallicity needing to be maintained. Vice-versa, for the low- 
mass galaxies, the metallicity stays low if the system is young and 
has a high amount of HI and SSFR. 

As observed in Fig.[T]and|2] within the range of stellar masses 
present at a given metallicity, the lower-mass end will still have 
higher SSFRs than the higher-mass end, but the difference com- 
pared to the lowest metallicities will be reduced. This is because 
the higher mass galaxies can have progressively more star forma- 
tion and still stay at mid-range to higher metallicities, while the 
lower-mass galaxies could have had a more bursty star formation 
in the past due to a higher amount of HI, which today would result 
in a high metallicity, low HI content and hence low SSFR. 

We can argue that more star formation increases the metallic- 
ity, which is sensible for massive galaxies, in which their SSFR in- 
creases as metallicity increases. It is noteworthy however, that this 



does not apply for low-mass galaxies, which follow the opposite 
behavior. 

Therefore, we propose that galaxies will follow a downsizing 
scenario in which the HI content is an important driver of their evo- 
lution. Low-mass galaxies with a large amount of HI will process 
their gas on longer timescales and thus show low metallicities to- 
day. On the other hand, low-mass galaxies with a lower HI might 
have evolved differently. These galaxies might have experienced a 
more bursty SF in the past that exhausted their gas and increased 
their metallicity. 

Aditionally, a combination of the infall of pristine gas and en- 
vironment could play an important role here. In this second sce- 
nario, infall can be more significant for low mass galaxies, dilut- 
ing their metallicity, and provide pristine fuel for high SSFRs. It is 
possible, even perhaps likely, that a combination of such scenarios 
are at play. To explain the observations presented here, and further 
explore the cartoon model proposed, more detailed simulations tar- 
getted at this problem are necessary. 
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